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We rigorously analyze the frequency response functions and antenna sensitivity patterns of three 
types of interferometric detectors to scalar mode of gravitational waves which is predicted to exist 
■ in the scalar-tensor theory of gravity. By a straightforward treatment, we show that the antenna 

sensitivity pattern of the simple Michelson interferometric detector depends strongly on the wave 
length Asgw of the scalar mode of gravitational waves if Asgw is comparable to the arm length 
of the interferometric detector. For the Delay-Line and Fabry-Perot interferometric detectors with 
arm length much shorter than Asgw, however, the antenna sensitivity patterns depend weakly on 
Asgw even though Asgw is comparable to the effective path length of those interferometers. This 
agrees with the result obtained by Maggiore and Nicolis. 
PACS numbers: 04.30.Nk, 04.80. Nn, 04.50.+h 
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I. INTRODUCTION 

To detect gravitational waves in the frequency band between ~ 100 and ~ 1000Hz, the TAMA300 has been 
constructed [jjj and TAMA300 team has already carried out several data taking runs ||. Other kilometer-size inter- 
ferometric detectors such as LIGO, VIRGO and GEO600 are also now under construction and will be in operation 
within a couple of years In addition, the Laser Interferometer Space Antenna (LISA) is now planned to be 

constructed in the next decade [|| to detect gravitational waves in the frequency band between ~ 10~ 4 and ~ 0.1Hz. 
Main purposes of these projects are (i) the detection of gravitational waves and (ii) the observation of relativistic 
astrophysical phenomena in a strong gravitational field. The purpose (i) is related to a fundamental question on 
the correctness of general relativity because the detection of gravitational waves could give a strong constraint about 
number of modes and propagation speed of gravitational waves . 

As an alternative candidate of general relativity, the scalar-tensor theory of gravity ||[)| has recently received a 
renewed interest. One of the reasons for its revival is that a theoretical approach to unification of interactions in terms 
of superstring theories suggests that the theory of gravity should be the dilatonic one (a kind of the scalar-tensor 
theor y) r ather than general relativity, i.e., the pure tensor theory |Tc| ]. 

In (Ly] , Shibata, Nakao and Nakamura investigated the gravitational collapse of a spherically symmetric dust ball 
to be a black hole in the framework of the Brans-Dicke theory (see also |l^] and jn|). They found that the amplitude 
of scalar gravitational waves can be 



or equivalently 



with the frequency ~ 100(10M©/M) = 10~ 3 (10 6 M©/M)Hz which is in an appropriate frequency range for kilometer- 
size detectors and LISA, respectively. Here, M, M Q , r, and oj are the mass of the collapsed object, solar mass, 
distance from the source to the detectors, and the so-called Brans-Dicke parameter, respectively. The noise level 
of the advanced kilometer-size detectors and LISA will be ~ 10 -23 for the most sensitive frequency band. Thus, 
although the present experimental limit is stringent ui > 3300 fl4|| , the amplitude of scalar gravitational waves may 
be high enough to be detected by the interferometric detectors if stellar mass black holes are formed in our galaxy or 
if supermassive black holes are formed in the cosmological scale. 

As pointed out in (T^| and reanalyzed by Maggiore and Nicolis |l5| , the antenna sensitivity pattern of interferometric 
detectors to the scalar mode is different form that of the tensor mode. This implies that if four interferometric detectors 
which are coincidently in operation detects gravitational waves, it is possible to ask whether the scalar mode exists or 
not using the difference of the arrival time and the difference of the antenna pattern among four detectors Jig] , Also, 
the peak frequency of the scalar mode is likely different from that of the tensor mode , suggesting that the two 

modes could be distinguished in the frequency domain even in the case of one interferometer. Thus, for exploration of 
the correctness of general relativity using interferometric detectors, it is important to clarify the antenna sensitivity 
pattern and frequency response function not only for the tensor mode but for the scalar mode of grav itational waves. 
From this motivation, analyses have been carried out in |TT1 fl5f| so far, but they were too crude [ fT!) or approximate 
p5[ . In this article, we carry out a rigorous reanalysis of the frequency response function and antenna sensitivity 
pattern of the interferometric detector to the scalar mode without any approximation: We directly solve the Maxwell 
equation for the propagation of a laser ray and the geodesic equation for motion of mirrors and a beam splitter in 
contrast to a previous work jlf|. Our straightforward treatment enables us to investigate the dependence of the 
antenna sensitivity pattern on the wave length of the scalar mode of gravitational waves without any ambiguity. 

This article is organized as follows. In Sec. II, we briefly review the linearized equation for gravitational waves in 
the framework of the scalar-tensor theory of gravity. In Sec. Ill, we solve the Maxwell field equation in the spacetime 
with the scalar mode of gravitational waves to derive a plane wave solution which describes the path of a laser ray 
in the interferometric detector. In Sec. IV, we present explicit solutions for timelike geodesies in the presence of 
scalar gravitational waves, since the mirrors and beam splitter in the interferometric detector move along timelike 
geodesies. We present the dependence of the frequency response function and antenna sensitivity pattern of Michelson 
interferometric detectors on the frequency of the scalar mode of gravitational waves in Sec. V, and of Delay-Line 
(DL) and Fabry-Perot (FP) interferometric detectors in Sec. VI. Sec. VII is devoted to a summary. 

We adopt the unit in which the speed of light is unity. We follow the convention of the metric and Riemann tensor 
in ]l9j . Throughout this article, we adopt the abstract index notation; The Latin indices a, denote a type of 
tensor, Greek indices components of a tensor, and two Latin indices, i and j, the spatial components of a tensor. 



II. GRAVITATIONAL WAVES IN THE SCALAR-TENSOR THEORY 

The field equation for the tensor part in the scalar-tensor theory of gravity is written in the form 



r _8tt , 



(V Q $)V fc $ - i«? ab (V c $)V c $] + -j-(V a V 6 $ - 5ah V c V c $), 



(2.1) 



where G a b and V a are the Einstein tensor and covariant derivative with respect to g a b, respectively. The function, 
w($), is the coupling function which determines the coupling strength between the gravitational scalar field, <I>, and 
the spacetime geometry. T a & is the energy-momentum tensor of matter or radiation. The field equations for the scalar 
field $ and matter are given by 



V C V C $ = 



1 



2lj 



du> 



(V C $)V C $ 



(2.2) 



and 



V a T ab = 0. 



(2.3) 
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In the Brans-Dicke theory, w is a finite constant and in general relativity, to = oo. Note that measurement of a 
gravitational lensing angle of a light ray from a distant quasar by the gravity of the Sun now constrains to be 
> 3300 (14|. 

Since the interferometric detectors are located in the wave zone of gravitational waves, we need to consider the 
linearized, vacuum (T ao = 0) field equations. We write the metric tensor and gravitational scalar field in the forms 



9ab — Vab + h a b, 

$ = $o(l-5$), 



(2.4) 
(2-5) 



where rj a b is the metric tensor of the background Minkowski spacetime and $0 is the value of <E> in the infinity. We 
assume that magnitude of each component of the perturbations h a b and (5$ is much smaller than unity. The linearized 
equations for the metric tensor and gravitational scalar field are written as 



^d c d c h ab - ^r) ab d c d d h cd + d c d {a h L b) = -d a d b 5<S> + r] ab d c d c 5$, 



and 



where 



d c d c 5<i> = 0, 



n b — 'I a cbi 
h a b = Kb - ^VabK: 

d c = r) cd d d . 

Introducing a new symmetric tensor field h a b defined by 

hab = Kb + VabS$, 



Eq. (2.6) is rewritten to 



- \d c d~h a b - \rjabd c ddh cd + d c d (a h c b) = 



(2.6) 

(2.7) 

(2.8) 
(2.9) 
(2.10) 

(2.11) 
(2.12) 



Equation (2.12) is completely the same as the linearized Einstein equation for the trace-reversed metric perturbations 
Jl9[ . This enables us to impose the transverse-traceless gauge condition on h a b as in general relativity as 



d a K = 0, 

Ko = 0. 

Under this gauge condition, the line element in the wave zone is written in the following form, 

ds 2 = (1 + 5$) {-dt 2 + (5 i:j + /^ T )*cW} , 
where hj? = hij is the transverse-traceless spatial tensor which satisfies 

d c d c hJ7 = 0. 



(2.13) 
(2.14) 

(2.15) 
(2.16) 



From Eq. (2.15), hj^ and 5$> are regarded as the tensor and scalar modes of gravitational waves, respectively, which 
constitute three independent modes of gravitational waves in the scalar-tensor theory. Hereafter, we refer to the scalar 
mode as scalar gravitational waves (SGW). 

In the rest of this article, we assume that the amplitude of the tensor mode hj? is much smaller than that of 
SGW for simplicity. Such a situation is realized when the source of the gravitational radiation is almost spherically 
symmetric. We also assume 



5$ = 5$(t + z), 



(2.17) 



in the coordinate system ( [2.15 ). 
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III. MAXWELL FIELD 



A trajectory of a laser ray is determined by solving the source-free Maxwell equation as 

g ac V a F cb = 0, 

where 

F ab = V a A b - W b A a = d a A b - d b A a . 



(3.1) 



(3.2) 



From Eq. (2.15), we find that the line element t akes conformally flat form in the coordinate system adopted here and 

(3.3) 



in the absence of the tensor modes. Then, Eq. (3.1) is rewritten to 

g ac V a F cb = (1 + ,5$)- V c d Q F cfc = (1 + 6$)- 2 (d a d a A b - d a d b A a ) = 0. 
Adopting the Lorentz gauge condition d a A a = 0, we obtain the equation for the vector potential A a as 

O a d a A b = 0. 



(3.4) 

Note that Eq. ( |3.4| ) is of the same form as the equation in the Minkowski spacetime. 

In the geometrical approximation, a laser ray is regarded as a monochromatic plane wave in the absence of SGW, 
i.e., 



A„ 



exp {ifl^t - c L • f)} , 



(3.5) 



where A^, Ol, c*l and f are the constant amplitude, constant angular frequency of the laser ray for the r est o bserver, 
constant direction cosine of the laser ray and position vector in the conformal 3D space, respectively. Eq. (3.4) im plie s 
that even in the presence of SGW, the components of the vector potential A M is completely the same as in Eq. ([3.5]), 
and the trajectory of the laser ray is straight in the conformal 3D space, i.e., 



dr 



(3.6) 



From the above result, one might think that the laser interferometric detector could not detect SGW since SGW do 
not affect the Maxwell field described in the conformally flat coordinate system. However this is the misunderstanding. 
We should take into account that the laser ray is repeatedly reflected by mirrors which are perturbed in the presence 
of SGW. This non-trivial perturbed motion causes variations of the vector potential A^ and consequently generates 
side-band waves which are just observed by interferometric detectors. 



IV. TIMELIKE GEODESICS 



The mirrors and beam splitter are suspended by wires so that they can freely move at frequencies much higher 
than the pendulum frequencies. For simplicity, we assume that the emitter of laser rays is also a free mass. Then, the 
t rajec tories of these parts of detectors are determined by solving the timelike geodesies in the spacetime described by 

Introducing the retarded and advanced time coordinates u and v, 

u = t — z and v = t + z, (4-1) 
equations for a timelike geodesic are given by 

r d + ^)?l=0, (4.2) 



dr y dr 

f'l + J*)^ } =0, (4.3) 



dr \ dr 
d f/-. ^.-.dv 



' M>)^ \ =0, (4.4) 



and 
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Hereafter, we assume that the mirrors and beam splitter are at re st i n t he a bsence of SGW. In this case, we can 
set 5& = using freedom to scale the coordinates. Integrating Eqs. (4.2)— (4.5), we obtain 



and 



dx dy 
dr dr 



(4.6) 



(4.7) 



where C v is an integration constant. From the normalization of the geodesic tangent vector and the above equation, 
we derive 



du 1 
dr C v 

Since dz/dr = and dt/dr = 1 for (5$ = 0, the integration constant C v is equal to 1. As a result, we find 

|-i-i»w, 



dz 
rf7 



1 



By integrating Eqs. fl4.6| )-( 4.10 ), the trajectory of a free mass is derived as 



V = Uu 
z = — - 



5$(v)dT = Zi 



t+Zi 



5&(v)dv, 



(4.8) 

(4.9) 
(4.10) 



(4.11) 
(4.12) 

(4.13) 



where x\, j/j and z\ are integration constants which specify the location of the free mass in this coordinate system 
(2.15) for t — > — oo, and we have used 



dv 



(4.14) 



Then, we derive a trajectory of a free mass for an arbitrary incoming direction of SGW, which is obtained by a 
spatial rotation of the coordinate system ( [2.15 ) as 







: 


-1 








cos v cos ip — sin ip sin v cos ip 
cos 9 sin ip cos ip sin 6 sin ip 
— sin 8 cos 8 



where the direction cosine c*sgw of SGW in the new coordinate system is 

csgw = ( — sin#cos(y3, — sin 8 sin ip, — cos 8). 
The explicit form of the trajectory of the free mass in the new coordinate system is written as 



x = x 



1 f 

i — - sin 8 cos ip / d$>(v)dv, 

2 J-oo 

1 f t+A 
i — — sin 8 sin ip I 6$>(v)dv, 

^ J -co 



y = y\ 



Z = Zi COS 

2 



t+A 



S$(v)dv, 



(4.15) 

(4.16) 

(4.17) 
(4.18) 
(4.19) 



where x\, yi and z\ are constants which specify the location of the free mass in this new coordinate system in the 
infinite past t — * — oo, and 

A = Xi sin 9 cos ip + y\ sin 9 sin ip + Z\ cos 9. (4.20) 

Since the spacetime metric takes conformally flat form in the absence of the tensor modes, we might think that 
SGW would be longitudinal wave. However, from careful investigation of the geodesic deviation equation, we find 
that the tidal force associated with SGW acts along the directions orthogonal to the propagation direction of SGW 
. From this point of view, SGW should be regarded as the transverse wave. It is also worthy to note that even if 
free masses are at rest in the present coordinate, the relative distances between them are not necessarily constant if 
SGW exist. Converse is also true. 



V. DETECTION OF SGW BY A SIMPLE MICHELSON INTERFEROMETER 



In the following, we identify the x and y axes with the directions of arms of the interferometric detector (cf, Fig. 
|l|). We hereafter refer to the arm on the x (y) axis as Arml (Arm2). 

In the Michelson interferometric detectors, a laser ray first reaches the beam splitter and then is divided into two 
rays; one enters into the Arml and the other into the Arm2. Each ray is reflected by the mirrors and then returns 
to the beam splitter again. The path difference of the two arms is chosen in such a way that the interfering light is 
dark at the anti-symmetric port for achieving the best signal-to-noisc ratio. Here, wc focus on two rays which finally 
reach x = on the beam splitter coincidently at t — (cf , Fig. |l|) . 

To calculate the side-band waves of laser rays, we need the trajectories of the beam splitter and mirrors only in the 



x-y plane. Using Eqs. (4.17)— (4.19), we obtain the trajectory of the center of the beam splitter as 



x — xb(£) = — sin#cos<^ / S$(v)dv, (5- 1 ) 



— oo 

t 



y = y-B (t) = - - sin 9 sin ip I 6&(v)dv. (5.2) 



2 
1 
2 

The trajectory of the mirror in Arml is given by 

1 f t+x 

x = XM(t) — L\ — - sinflcos^ / 5&(v)dv, (5.3) 

* J-oc 

where 

X = L\ sin cos 93. (5-4) 

Note that we do not need the y-coordinate of this mirror for calculating the side-band waves since the motion in the 
y-direction does not affect the path difference in the linear order in 5$. From essentially the same procedure, the 
trajectory of the mirror in Arm2 is given by 



1 



t+Y 



y = yu(t) = L 2 - -sin 9 sin p 5$(v)dv, (5.5) 
^ J -00 

where 

Y = Li sin^sin^. (5-6) 

In this case, it is not necessary to know the x-coordinate of this mirror. 
The trajectory of the emitter is written as 



1 



t+x E 



x — xn(t) — — £e — — sin#cos<y9 / 5&(v)dv, (5.7) 
where 

Xe = — -^e sin 6 cos ip. (5.8) 
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Information of the y-coordinate of the emitter is also not necessary. 

Since we regard the emitter as a free mass, we need to carefully treat the phase of the laser ray. The amplitude A 
of the laser ray at the emitter is written as 

A = Ae ,sllTE , (5.9) 

where Ah and are constant and te is the proper time of the emitter. The relation between the coordinate time t 
and proper time te is given by Eq. (4.9). Integrating this equation, we obtain 

i rt+X E 

T E (t) =t+ - / 5$(v)dv. (5.10) 
Thus, the amplitude A is expressed as a function of t as 

rt+X E \ ~| 

A(t) = Ah exp I iU L I t + ~ I S$(v)dv } . (5.11) 




Then, we consider the ray entering into Arml. Tracking from t — t B toward the past, the laser ray reached the 
mirr or in Arml at t = t-g, — L\ + 0(8$), and then it reached the emitter at t = t B — 1L\ — £e + 0(5$>). As shown in 
Eq. (3.5), the plane wave solution of the Maxwell equation is completely the same as that of the Minkowski spacetime 
in the conformally flat coordinate system. This implies that the amplitude A\ of the laser ray at the beam splitter 
agrees with the amplitude A at the emitter at t = t s — ri , where r% is the coordinate path length of this ray from the 
emitter to the final beam splitter. Using Eq. ( 5.1l| ), we obtain the amplitude A\ of the laser ray at the beam splitter 
on the output port side at t = as 

in L {t B -r 1 + -j 5®(v)dv\ y (5.12) 

where p& is a reflection coefficient from the detector side (— from the laser side) and Uq is transmission coefficient 
of the beam splitter. 

From the mirror at t — £b — L\ + 0(<5$) to the beam splitter at t = the coordinate path length of the laser ray is 
equal to xm^b — Li) + 0(5& 2 ), and from the beam splitter at t = te — 2Li + 0(8$) to the mirror at t = t-Q — Li + 0(5$), 
the coordinate path length is xm(^b — L\) + 0(5$ 2 ). From the emitter at t = t-Q — 2Li — £b + 0(5$>) to the beam 
splitter at t — t-g — 2Li + 0(6$), the coordinate path length is — xe(^b — 2L\ — £e) + 0(S$ 2 ). Thus, n is given by 



ri = x M (tB - Li) + x M (tB - Li) - x E (t B ~ 2Li - £ E ) 

rtB-Li+X / .t B -2Li-^E+^E S 




6<f>(v)dv. (5.13) 



Substituting the above expression into Eq. ( 5.12 ), the amplitude A\ is written in the form 

A x = p B ^Ae ,!!l(iB_2Il ^ E) {1 + AAx + 0(5t> 2 )} , (5.14) 

where 

AAi = I sin 9 cos ( 2 y -J \+J \5<P(v)dv. (5.15) 

From the same procedure as that for the laser ray entering into Arml, the amplitude Ai at the beam splitter on 
the output port side for the laser ray entering into Arm2 is expressed in the form 

itt L lt B -r 2 + -J 5$(v)dv\\, (5.16) 

where r-i is the coordinate path length of this ray from the emitter to the final beam splitter. 

As shown in Fig. [I], the trajectory of the laser ray considered here is restricted on the straight line x = when it 
goes through the Arm2. Thus, this laser ray comes back to the beam splitter first at 
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y = y B (t B -2L 2 )-x B (t B -2L 2 ), (5.17) 

and finally at 

V = 2/b(*b) - zb(*b)- (5.18) 

From the same procedure as in the case of Arml, the coordinate path length r 2 from the emitter to the final beam 
splitter is given by 

r 2 = 2/m(*b - L 2 ) - {jte(is) - x B (ts)} 

+ 2/m(*b - ^2) - {j/B (^b - 2i 2 ) - x B (t B ~ 2L 2 )} 
-x E (t B -2L 2 -£ E ) 

ftB rt^—2L2 ftB — 2L 2 — £e+^E 

2i 2 + — - sin ^ cos <p 



—00 «/ —00 



sin<^ [ / + / }6${v)dv. (5.19) 



t B -2L 2 



Substituting the above equation into Eq. (5.16), we obtain 

A 2 = -p B VBA h e m ^ tB - 2L2 -^ {1 + AA 2 + 0{S^ 2 )} , (5.20) 
at the final beam splitter, where 

• J ( ftB ftB—2L 2 ptB—^Li—tB+X-B 

AA 2 = — Ol < sin # cos <^ 



00 J — 00 



+ sin^sin^ ( / +/ +/ U$(u)du. (5.21) 

\Jt B JtB-2L 2 J J-00 J 

To maintain the output port to be dark in the absence of SGW = 0, the following condition should be satisfied 

fl L (L 1 -L 2 )=mn, (5.22) 
where m is integer. Therefore, in the presence of SGW, the side-band waves are generated with the amplitude 

A out = Ai + A 2 

= p B i^ B A L e in ^- 2L ^\AA 1 - AA 2 ) 

r+00 

= PB u B A h n h e m ^- 2L ^ 5$(u;)e iutB H M (u)dLJ, (5.23) 



where 



1 



+ OC 



8$(u) = — / e- wv 5$(v)dv, (5.24) 

and Hyi{uj) is a frequency response function to SGW of Michelson interferometric detectors. (Note that hereafter ui 
does not denote the Brans-Dicke parameter but the angular frequency.) It is worthy to note that difference between 
the coordinate time t B and the proper time at the detector is 0(<5$) so that it only contributes to a higher order 
correction for the side-band waves. 
The explicit form of Hm(oj) is 

H m (lj) = ^-sin6»|cos^(2e i " il ( sinecosv - 1 ) - 1 -e" 2iwil ) 

-sin^(2e ML2(sinesin¥ '- 1) -l-e" 2 ^ 2 ) \. (5.25) 
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For simplicity, we restrict our attention to the case 



L 1= L 2 = L. (5.26) 
Then, in the limit uj — > 0, the response function becomes 

lim H M (uj) =iL sin 2 9 cos 2ip, (5-27) 



UJ^O 



which agrees with the antenna sensitivity pattern obtained by Maggiore and Nicolis [ ^5[ . However, for SGW of general 
frequency uj, the antenna sensitivity pattern is different from Eq. (5.27) (cf, Fig. ||). 



To see the dependence of the amplitude of Hm(uj) on the frequency uj of SGW, we pay particular attention to the 
incident direction 9 = n/2 and ip = 0. Then, we obtain 

|HmM| 2 = ^^. (5.28) 

Thus, Michelson interferometric detectors are sensitive to SGW of uj < u>m where 

uj m = n/(2L). (5.29) 



From Eq. (5.27), we find that the longer the arm length L is, the better the sensitivity of Michelson detectors is. 
However, longer arm length leads to smaller upper bound ujm of the sensitive frequency range. To optimize the laser 
interferometric detector against SGW of an aimed frequency, the arm length of the laser ray should be equal to a 
quarter of the wave length of SGW. 

The sensitivity patterns on 9 = n/2 plane for the frequency uj = ujm of aimed SGW are shown in Fig. ^(a). 



Obviously, the sensitivity patterns for uj = ujm and uj = wm/2 are different from Eq. (5.27). On the other hand, for 
u j — o ;m/100; *.e-, when the wave length of SGW is 400L, the sensitivity pattern is almost the same as that in Eq. 
fl5~27l ). These consequences imply that the sensitivity pattern of Michelson interferometric detectors depends strongly 



on the frequency in the neighborhood of uj = ujm- 

VI. DETECTION OF SGW BY DELAY-LINE AND FABRY-PEROT INTERFEROMETERS 

As explained in the previous section, the arm length of Michelson interferometric detectors should be equal to a 
quarter of the wave length of SGW to optimize the sensitivity. The wave length of SGW generated by a spherical 



gravitational collapse to a black hole, however, could be ~ 100(M/M Q )km for stellar mass black holes |11|. It is in 
practice difficult to construct a laser interferometric detector of such long arm length on the ground. Due to this 
reason, DL or FP interferometric detectors in which it is possible to make the effective optical path length much 
longer than the arm length have been adopted as ground-based detectors in practical projects. In this section, we 
study the frequency response function and antenna sensitivity pattern of DL and FP laser interferometric detectors. 

Following the notation in Sec. V, we identify the x and y axes with the direction of the arms of the interferometric 
detector (cf, Fig. [|). We again refer to the arm on the x (y) axis as Arml (Arm2). As in the case of Michelson 
detectors, we need the trajectories of t he m irrors and beam splitter only in the x-y plane. The trajectory of the center 



of the beam splitter is given by Eqs. ( J5.l[ ) and (5.2). The trajectories of the front mirror and of the end mirror in 
Arml are given by 

x = XMf(t) = t\ — — sintfcost^ / 6§(v)dv, (6.1) 

2 J-oo 



1 



t+X c 



where 



x — XMc(t) = L\ + £i — — sm9 cost^ / 5$(v)dv, (6.2) 



X{ = l\ sin 9 cos ip, (6.3) 
A = (Li +£i)sin6»cos^. (6.4) 



The trajectories of the front mirror and the end mirror in Arm2 are given by 
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y = VMf(t) =£ 2 --^sux9 sin <p 



t+Yf 



5$(v)dv, 



y = VM e (t) = L 2 + £-2 - 2 sin ^ sin 



t+y c 



5$(v)dv, 



where 



Yf = £2 sm9 shxip, 

Y c = (L 2 + £2) sin 9 sin ip. 



(6.5) 
(6.6) 



(6.7) 
(6.8) 



We do not need to calculate the y-coordinates of mirrors in Arm 1 and the ^-coordinates of mirrors in Arm2 to obtain 
side-band waves since the motion of these mirrors in these directions affects the path difference from the seco nd order 
in <5$. The t rajectory of the emitter and the amplitude A of the laser ray at the emitter are given by Eq. (5.7) and 
by Eq. ( 5.11 ), respectively, as in Michelson detectors. However, the definitions of L% and L 2 are different from those 
in Michelson detectors (cf, Fig. 0). 



A. Delay-Line Interferometer 

As in Michelson interferometric detectors, a laser ray in DL detectors first reaches the beam splitter and then is 
divided into two rays; one enters into the Arml and the other into the Arm2. In the DL detectors, however, the 
ray entering into the regions between the front and end mirrors is reflected by these mirrors by (2N — 1) times until 
returning to the beam splitter again to increase an optical path length by 2N times as long as the arm length. This 
mechanism changes the response of the detector. 

As in Sec. IV, we pay attention to two rays which finally reach the beam splitter coincidently at t = t& (cf, Fig. 
|3j). First, we consider the ray entering into Arml. For notational convenience, we introduce 

*f(n) =t B -{4 + (2n-2)Li}, (6.9) 
t e (n)=*B-{4 + (2n-l)L 1 }, (6.10) 
*E(n) = *b - (24 + 4: + 2nL t ) . (6.11) 

From the same procedure as in Michelson interferometric detectors, the amplitude A\ of the A^-fold ray at the beam 
splitter on the output port side is expressed in the form 



1 



A x = pb^b^l exp liCl h \ t B - r (N) + -J S1>(v)dv j j , (6.12) 

where the definition of pb and v-q are the same as for the Michelson detector, and rvjv) is the coordinate path length 
of this ray from the emitter to the final beam splitter. 

In the following calculation for the coordinate path length 77jv)j we assume that the trajectory of the ray is parallel 
to the x-axis. Strictly speaking, the ray is not parallel to the a;-axis in a real DL interferometric detector, but the 
effect of it on the path length is the second order in <5$. Then we obtain 

JV 

r (N) = 2>Mf(*f(l)) + ^ {2x Mc {t c (k)) - x Mf{h(k)) ~ ^Mf(*f(fc+1))} + ^Mf(*f(JV+l)) _ £e(£e(AT)) 
fc=l 

= 2£ 1 +£ E + 2NL 1 + AR {N) , (6.13) 

where 

AR {N) = --sm9costpl + - +J2[ +/ >6$(v)dv. (6.14) 

2 [J-00 J-oa J-00 k=1 \ Jt m +X t Jt t(k + 1) +X { J J 

Substituting Eq. (6.13) into Eq. ( 6.12| ), the amplitude A\ becomes 

A x = PB is B A h e tn ^ tB - 2ei - eE - 2NL ^ {l + AAi + 0(<5$ 2 )} , (6.15) 

where 
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(1 ftE(N)+XE \ 

AR {N) --J 5<5>(v)dv\. (6.16) 

Next we consider the ray entering into Arm2. As above, we introduce 

i f{n) =t B -{e 2 + (2n-2)L 2 }, (6.17) 

*e(n) = *b - {4 + (2n - 1)L 2 } , (6.18) 

*B(n) = *b - (2^2 + 2nL 2 ), (6.19) 

*E(n) = *b - (2£ 2 +Ie + 2nL 2 ), (6.20) 

and in addition, 

w B {t) = y B {t) - x B {t). (6.21) 
The path length f(jv) of the iV-fold ray from the emitter to the final beam splitter is given by 

AT 



>6®(v)dv 



r(N) = J/Mf(if(l)) - W B (t B ) + ^ { 2 2/Mc(^o(fe)) - J/Mf(*f(fc)) - J/Mf(if(fe+1))} + J/Mf(if(iV+l)) - Wb(*B(JV)) - Ze(*E(JV)) 

k=l 

= 2£ 2 +£ E + 2NL 2 + AR (N) , (6.22) 

where 

1 f /•*f(N+i)+>f ri f (i)+Yi N / ,f e(fc) +r e rt c(k) +Y a 

AJ? (JV) = --sintfsin^ / + +^ / +/_ 

(/•*E(1V)+^E /-*B / , *B(iV)\ 

/ - - }6$(v)dv. (6.23) 

From the same procedure as in the case of Arml, the amplitude A 2 of the laser ray entering into Arm2 is 

^2 = -p B ^B^Le inL( * B - % - fe - 2JVi2) {1 + AA 2 + 0(<5$ 2 )} , (6.24) 
at the final beam splitter on the output port side, where 

^ / , *e(n)+-X'e 



1 / / , * E ( W )~'~-^ E /"* B 

+ - sin cos ( ' 



AA 2 = -ifl L ^AR {N) - - J 6<S>(v)dvj . (6.25) 

To maintain the output port to be dark in the absence of SGW, the following condition should be satisfied 

nj J (l 1 +NL 1 )=n h (l 2 + NL 2 ) + miT, (6.26) 

where m is integer. Thus, in the presence of SGW, side-band waves are generated with the amplitude 

Aut = Ax + A 2 

= p B ^ B A L e m ^- 2t ^- 2NL ^(AA 1 - AA 2 ) 

/+oo 
S^>(u;)e luJtB H Bh (u;)du;, (6.27) 
-oo 

where -Hdl(w) is a frequency response function to SGW of DL interferometric detector. The explicit form of J?dl(w) 
is 

1 



sin0cos(^je^ Xf - £l - 2ArLl ) _ i + e M^-<0 _ e -«o,(/ a +iVL a ) 

C E -£ E ) ( e -2ioj(i 2 +NL 2 ) _ e -2iLo(£ 1 +NL 1 )^ 

-liuLx l-\ _ -2injNL 1 \ } 

, £ V 1 e )_,2 e iu{X a -li+L 1 ) _ e iu(Xt-l 1 +2L 1 ) _ ia)(X f -<i)\ I 

J _ g-2iwLi f 
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sin sin </? < e 



,iu(Yt-t 3 -2NL 3 ) _ 1 + e iu)(Y f -e 2 ) _ e -2iu>{l 2 +N L 2 ) 



—2iu>L2f-\ _ „—2iuiNL 3 \ 

Z V 1 e l(2 e iu(Y -e 2 +L2) _ e iu(Y{-l 2 +2L 2 ) _ e iu{Y { -t2)\ 

\ _ e -2iuL 2 <> ' 



(6.28) 



For simplicity, we restrict our attention to the case 

£ 1 =l 2 =£ and Li = L 2 = L. (6.29) 

In the practical DL interferometric detectors, I <C L and loL < 1 for the most sensitive frequency band of the detector. 
Then, for u>£ <C 1, the frequency response function of DL detector becomes 

HuL (u) ~ sin (^) e -to(iy-i)L sin J cos y ( 2e ^(-fa» « c,s v -x) _ 1 _ e -2i uL) 
2wsin(wL) l_ 

- sinv3(2e ML(sinf,sinv " 1) - 1 - e -^ L ) 

= ^NI^ e _ iu{N _ x)L 
sin(wL) 

where Hyi{ui\ L) is the frequency response function of Michelson interferometric detector of the arm length L\ — Li — 
L. Hence, the antenna sensitivity pattern of DL interferometric detectors is the same as that for Michelson detectors 
of the arm length L, except for the overall factor depending on u>. In the limit lu — > 0, the frequency response function 
for two types of detectors coincides as 

lim # DL (w) = iNL sin 2 6 cos 2ip. (6.31) 

However, for finite lu, this is not the case. To see the dependence of the amplitude of Z/dl 011 the frequency to of 
SGW, we pay particular attention to the case that the incident direction of SGW is — ir/2 and ip = 0. In this case, 
we find 

This implies that the sensitivity of DL interferometric detector of iV-fold is high for oj < wdl, where 

wdl = n/{2NL). (6.33) 

To optimize DL interferometric detectors against SGW of an aimed wave length AsgWj the folding number N is usually 
set to be Asgw/4Z so that cjdl = 27t/Asgw is satisfied. As mentioned above, in the case of practical interferometric 
detectors constructed on the ground, N is much larger than unity. Since the dependence of the antenna sensitivity 
pattern on the frequency of SGW is the same as that of Michelson one, and since the DL detector of N 3> 1 is 
sensitive for the wave length Asgw > 27r/tJDL 3> L, the antenna sensitivity pattern of DL for aimed SGW is almost 
the same as Eq. ( |6.31 ), even if the wave length of SGW is comparable to the optical path length. In other words, we 



can consider that the amplitude of SGW is effectively constant in computing the response of the detectors. 



B. Fabry-Perot Interferometer 



The FP interferometric detectors form the Fabry-Perot cavity between the front and end mirrors. We define that 
the reflection coefficients of the front mirrors for the light ray incident from the FP cavity side and laser side are 
Pm and — pm, respectively, and that the transmission coefficients of these are vm- Hereafter, we consider a loss- free 
system in which 

"Af + /& = !• (6-34) 
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We assume that the end mirrors have a perfect reflectivity, i.e., the reflection coefficients of those are —1. On the 
other hand, the reflection and transmission coefficients of the beam splitter are assumed to be pe and v B as before. 

As in the case of Michelson and DL interferometric detectors, we focus on laser rays which finally reach x = on 
the beam splitter at t — ts- First, we consider the ray entering into Arml. The path length r^o) of the ray which 
does not enter into the FP cavity is 

r (0) = 2x M f(£f(i)) - £e(*e(o)) 

= 2h +£ E - ^ sin (9 cos ^ I 2 J -J I 5&(v)dv. (6.35) 

From the same procedure as in the case of DL interferometer, the path length r(„) of the ray which enters into the FP 
cavity and is reflected by the front and end mirrors by (2n + 1) times until returning to the beam splitter is given by 

r(n) = f(o) + 2nLi + Ar (n) , (6.36) 

where 



Ar (n) =-- sin0 cos <pl / -/ + £ ( / +/ )\5®(v)dv. (6.37) 

[•^ *f(l)+-ff •'*E(0)+^E fc = J \"'tf(fc)+Jff •> tt(h+\)+Xf J J 

As in the derivation of Eq. ( 6.12| ), the amplitude ^4i(o) of the path length r( ) at the beam splitter on the side of the 
output port is obtained as 

A m = -wPBPMe in ^- 2t ^- Sr ^\ (6.38) 

where 

^ | /"*f(l)+-Xf r*E(0)+^-"E | 

fr( ) = — <2sin0cos<W + (1 - sinflcos^) / > <5$(w)d-y. (6.39) 

2 I i-oo J-oo J 

The amplitude Ai( n ) of the path length r(„) at the beam splitter on the side of the output port is given by 

A 1{n) = -(^) A m {~ PM e- 2in ^) n e- in ^\ (6.40) 
\PmJ 



^ /•*E(n)+-^E 



where 

<5r (n) = Ar (n) - | / (6.41) 

•'*E(0)+-^E 

Hence when the ray reaches the beam splitter on the side of the output port at t = t B , the amplitude A\ of the laser 
ray is 

oo 

Ai = ]T A 1(n) = ~^ BPB p M A L e^ tB - 2e ^( ai + A Ax) + 0(5<S> 2 ), (6.42) 

n=Q 

where 

2 -2iO L L! 

at = l+ — j-M ^y-, (6.43) 

AAx = -ztt L |aifr (0) - (j^) fj(- PM e- 2lOLLl ) n fr (n) | . (6.44) 

Next, we consider the ray entering into Arm2. From the same procedure as that for the Arm 1, we obtain the path 
lengths f (o) and f („) as 

f( ) = 2j/Mf (tf(l)) - ^B(te) - wb (£b - 2^ 2 ) - £e(*e(o)) 

if/ /"*f(l)+^f /"*f(i)+^f\ / /"'b rt-B~2t 2 ri-E(o)+ x E.\ | 

= 2£ 2 +£ E --sinl^sm^ / +/ \+aaatp\l + \6${v)dv. (6.45) 

^ \Jt B Jt B -2e 2 I \J-00 J-00 J-oa I I 
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r(n) = r (0 ) + 2nL 2 + Ar (n) , (6.46) 

where 



Af (n) = 2 



-sinflsuW / + / + / + / >(5$(w)c?w 

2 [^((ij+Yf Jt B( „) fe=1 \"'if( fc )+V f Jt Kk+1) +Y{) J 



1 / ptB—2^2 / i *E(T2)+^E\ 

+ - sin 6> cos (p I / +/ J 5$(v)dv. (6.47) 

Here, the path length f(„) of the ray is computed for the ray which enters into the FP cavity and is reflected by the 
front and end mirrors by (2n — 1) times until returning to the beam splitter. 
The amplitude A 2 of the laser ray entering into Arm2 is written in the form 

A 2 = v BPBPM e m ^^-^\a 2 + AA 2 ) + 0(<5$ 2 ), (6.48) 

where we assume that the ray reaches the beam splitter on the side of the output port at t = ts, and 

2 -2iO L L 2 

AA 2 ee -itt L \a 2 Sf (0) - (^-) Y,{- PM e- 2 ^ L -) n 8f (n) ] , (6.50) 



' n— 1 ) 
tf(i)+*f /•*f(D+ y f\ / /-*B rt-z-lli 



- - < sin sin ip I / + / + sin # cos I / + / 

2 WtR Jtn-2£ 2 / W-oo J- 



£b Jtu— 2£2 J \J—oo J— oo / 

/*E(0) +-f E | 
U$(v)cfo, (6.51) 

<yf (n) ee Af (n) -\f_ () 5$(v)dv. (6.52) 

^*E(0)+^E 

We assume that the following resonance condition holds in the absence of SGW, 

2n L L, = (2n 4 + 1)tt, (6.53) 
where rii is integer. Since the output port should be dark in the absence of SGW, i\ and l 2 should satisfy 

Cl L (£ 2 =rmr, (6.54) 
where m is integer. Hence the output in the presence of SGW is given by 

A out = A 1 + A 2 = -^pBpuA^^-^-^iAA, - AA 2 ), 

/+oo 
S$(u;)e^ tB H F p(Lu)dLu, (6.55) 
-oo 

where Hpp(ui) is the frequency response function of the FP interferometric detector to SGW. 
For simplicity, we focus on the case 

l x = t 2 = i and L X =L 2 = L. (6.56) 

Then the response function Hpp(cu) is written as 

fr F p( w ) = --L 



sin<9cos^|(2a - b)e iu( - Xt ~ e) -(a- &)e" 2 ^ - a 

lb U e MXa-i) _ e iw(x f -e+L) _ e iv(Xi-e-L)\ 1 
2sin(wL) V ) J 



sin6»sin^|(2a - b)e ihj{ ~ Y ^ -{a- b)e 



2iul _ 



lb U e iu(Y -t) _ e ioj(Y{-t+L) _ e iu>(Y t -e-L) 

2sin(cjL) V 



(6.57) 
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where 



a = l- h M y (6-58) 
Pm(1 - Pm) 

6= UJ (w^-T^j- (6 ' 59) 

Due to the same reason as in the case of DL interferometric detectors, we assume that £ is much smaller than L. 
Then for lu£ <C 1, the frequency response function H-pp(u) is derived as 

ibe iuL 

H fp (lj) ~ - , r . H M (cj; L). (6.60) 

Thus, FP interferometric detectors also have the same angular dependence of the frequency response function as 
Michelson detectors of the arm length L. In the limit ui — > 0, the response function H F p(uj) becomes 

lim H FP (cu) = '^ + PM \ Lsin 2 6»cos2(^. (6.61) 
w^O p M (l - p M ) 

In the case of finite oj, the absolute value of the response function |7?fp(w)| for 9 = tt/2 and <p — becomes 



( 1 + p M ) / 1 - cos(2ojL) 

p M w V 2{l + p^ ~2p M cos(2cji)}' 



\H FP (uj)\ = — —\lTrn-, — 2 o„ ZZZTK-TTvT - ( 6 - 62 ) 



We show \Hpp\ for various pm in Fig. [|. Note that \Hpp \ depends on the value of the reflection coefficient pm of the 
front mirror sensitively. For example, in the case of TAMA300, pm is equal to 0.988. The interferometric detector 
such large pm is sensitive for to < 10~ 2 7r/(2L) for which the wave length is much longer than that of the arm length 
of the interferometer. Hence the antenna sensitivity pattern of such a FP interferometric detector is almost the same 
as Eq. ( |6.6l| ). This result is consistent with Maggiore and Nicolis [ p| . 



VII. SUMMARY 



We have rigorously analyzed the sensitivity of interferometric detectors to SGW. Since we calculate the side-band 
waves of laser rays directly solving the Maxwell field equation, we were able to obtain the dependence of the antenna 
sensitivity pattern on the frequency of SGW without any ambiguity. As a result, we found that the antenna sensitivity 
pattern of the simple Michelson interferometric detector depends strongly on the frequency of SGW. These features 
are essentially the same as in the case of tensor mode of gravitational waves ^(J . LISA is the interferometric detector 
and will be simply Michelson. Thus, the dependence of the antenna sensitivity pattern on the wave length of SGW 
has to be taken into account as in the case of the tensor modes. 

We have also found that the dependence of the antenna sensitivity patterns of both the Delay-Line and Fabry-Perot 
interferometers on the frequency of SGW is proportional to that of the simple Michelson one. This implies that if 
the wave length of SGW is comparable to the arm length of the interferometric detector, the antenna sensitivity 
patterns of Delay-Line- and Fabry-Perot interferometric detectors are different from the result obtained by Maggiore 
and Nicolis who assumed that the amplitude of SGW is constant in their calculation. However, the effective optical 
path length of Delay-Line- and Fabry-Perot detectors will be optimized for SGW of the wave length much longer 
than the arm length. In this case, the wave length of SGW is much longer than the arm length of the interferometric 
detector and the amplitude of SGW can be considered to be effectively constant for evaluating the response, even 
if the wave length of SGW is comparable to the optical length. Therefore, the antenna sensitivity pattern is almost 
same as the results obtained by Maggiore and Nicolis fll5||. 
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FIG. 1. The definitions of various quantities for Michelson detectors, (a) A configuration in the absence of gravitational 
waves, (b) A schematic configuration in the presence of scalar gravitational waves. 
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FIG. 2. (a) The frequency response functions of Michelson interferometric detectors for to = com, wm/2 and wm/100 on 
9 = 7r/2 as a function of (p (the solid, dashed and dotted lines, respectively), (b) The frequency response function of the same 
detector for the same w as (a) on ip — is depicted as a function of 6. In both cases, we take the units in which L\ and L2 are 
unity. 
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FIG. 3. The definitions of various quantities for DL and FP interferometric detectors, (a) A configuration in the absence of 
gravitational waves, (b) A schematic configuration in the presence of scalar gravitational waves. 
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FIG. 4. The frequency response function of FP interferometric detectors at 6 = n/2 and 95 = as a function of ui. In this 
figure, we adopt the units in which L is unity. The solid, dashed, and dotted lines denotes for pM = 0.988, 0.900 and 0.770, 
respectively. For pM = 0.770, the same sensitivity becomes almost same as that for Michelson interferometric detectors. 
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